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Based on the dielectric continuum model and transfer-matrix method, the completing polar optical phonon 
states in a wurtzite GaN-based superlattice (SLs) with arbitrary-layer complex bases are investigated. It is 
proved that 2™ types of phonon modes probably exist in a wurtzite nitride SL with n-layer complex bases. The 
analytical phonon states of these modes and their dispersive equations in the wurtzite GaN/AUGai-xN SL 
structures are obtained. Numerical calculations on a three-layer GaN/Alo.isGao.ssN/AIN complex bases SL 
are performed. Results reveal that there are interface optical (IO) phonon modes of one type only and four 
types of quasi-confined (QC) phonon modes in three-layer GaN/Alo.isGao.ssN/AIN complex bases SLs. The 
dispersive spectra of phonon modes in complex bases SLs extend to be a series of frequency bands. The 
behaviors of QC modes reducing to IO modes are observed. The present theoretical scheme and numerical 
results are quite useful for analyzing the dispersive spectra of completing phonon modes and their polaronic 
effect in wurtzite GaN-based SLs with complex bases. 
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1. Introduction 

Since the superlattices (SLs) notion was first brought forward nearly four decades ago [![, the 
electronics and opto-electronics properties of semiconductor SLs have attracted considerable at- 
tention due to their fancy miniband effect and a lot of excellent physical properties P, 0| • Among 
them, semiconductor SLs based on group-Ill nitride GaN, A1N, InN, and their ternary compounds 
Al^Gai-^N, In-rGai-zN have invoked special interest both in theoretical and experimental investi- 
gations [3|-l7|. This is mainly ascribed to the following evident characteristics: the nitride materials 
with widely and adjustable direct-band gap, strong atomic bonding and high carrier mobility make 
them attractive materials as a basis for the creation of reliable high-temperature and high-frequency 
electronic equipment and short- wavelength optoelectronics devices [5H7J. 

In order to further improve the performance of microelectronics devices based on SL structures, 
people brought forward the complex-bases SL, i.e., n layers per period (n > 2) [8HlO|. In contrast to 
the binary (two-layer per period) SLs, the complex bases SLs provide more degrees of freedom for 
the adjustable structural parameters, the miniband and minigap widths due to the introduction of 
additional layers in each SL period. This is very important for several important applications, such 
as infrared photodetectors, effective- mass filtering, and tuning of the tunneling current [Tll - ll3| . The 
previous research on the complex bases SL is mainly focused on the electronic structures and sur- 
face polariton properties in GaAs-based semiconductors [l^ - [l9| . But there are few investigations 
of the other physical properties, such as the crystal-lattice dynamical feature and optoelectronic 
characteristics in GaN-based SL with complex bases [2(| . Indeed, crystal-lattice vibrations (phonon 
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modes) have great effect on the electronic and impurity states, the excitonic recombination and car- 
rier transport, especially in low-dimensional quantum systems of polar semiconductors pH l22j|. At 
room and higher temperatures, the electron-phonon interactions and scattering by optical phonons 
play a dominant role effecting various properties of polar semiconductor quantum heterostructures, 
including hot-electron relaxation rates, interband transition rates, room-temperature exciton life- 
times, and many other optical and transport properties [23|,[24j|. Moreover, since nitride material 
usually crystallizes in the anisotropic hexagonal wurtzite structure, the crystal dynamical prop- 
erties of nitride SLs are more complex due to anisotropy of the crystal structure in contrast to 
those of cubic crystals (25l - l32j . Hence, it is necessary and important to investigate the polar optical 
phonon modes in wurtzite GaN-based SL with complex bases. 

On the basis of the dielectric continuum approximation and Loudon's uniaxial crystal model [33] , 
several authors have made their great contributions to the study of the polar optical phonons and 



their electron-phonon interactions in wurtzite nitride binary SLs [20, 341. 135l|. For example, Gleize 
and his coworkers [34j investigated the anisotropy effects on polar phonon modes in wurtzite 
GaN/AIN SLs, and the quasi-confined (QC) phonon modes and interface-optical (10) phonon 
modes were found in the SL structures. Medeiros, Anselmo and their cooperators studied the 
confinement of polar optical phonons in two-layer periodic GaN/AIN SLs [35| based on the dielec- 
tric continuum model (DCM), and the phonon dispersive spectra of the systems were displayed 
in their discussions. From the viewpoint of experimental investigations, Darakchieva et al. [36| 
studied the phonons in AIN/GaN SLs using a combination of infrared spectroscopic ellipsometry 
and Raman scattering spectroscopy. And the Ei(TO), Ai(LO) and E2 localized, and Ei(LO) and 
Ai(TO) delocalized superlattice modes were identified. Recently, Davydov et al. studied the 
lattice dynamics and Raman spectra of strained hexagonal wurtzite GaN/AIN and GaN/AlGaN 
SLs, and thus the 10 and the QC phonon modes were observed. Gleize and coauthors [38( probed 
the confined phonons in hexagonal GaN-AIN artificial SL structures based on non-resonant Raman 
scattering, and the QC optical phonon modes were observed. Gleize's group [39] also measured the 
angular dispersion of polar 10 and QC phonons in a hexagonal nitride SL, and the experimental 
data were in good agreement with the results of a previous calculation based on a DCM. Dutta 
et al. [40j reported Raman scattering results in binary SLs of GaN/(In)AlN, which also proves 
the availability of DCM for the description of polar optical phonon modes in wurtzite nitride het- 
erostructures. Therefore, the DCM and Loudon's uniaxial crystal model are reliable and will be 
adopted to the investigation of the optical phonon modes in wurtzite nitride SL with complex bases 
in the present paper. 

The main accomplishments and significance of this work lie in the following three points, (i) 
Via solving the Laplace equation of phonon potentials in wurtzite SLs with n-layer complex bases, 
it is confirmed that 2™ types of polar phonon modes, i.e. the propagating (PR) mode, the 10 mode, 
the (2™ — 2) QC modes probably coexist in wurtzite complex-base SL structures, (ii) By using the 
transfer- matrix method [1,0], the polar optical phonon states and corresponding dispersive equa- 
tions for the 2™ types of polar phonon modes in wurtzite nitride complex-bases SLs are given within 
the framework of the DCM and the Loudon's uniaxial lattice model, (iii) As an example, a wurtzite 
GaN/ALGai_a;N/Alj,Gai_j,N SL with three-layer complex bases is chosen to check the theoreti- 
cal scheme described in the present paper. The numerical result reveals that there are five types 
of phonon modes in the wurtzite GaN/AL r Gai_. E N/Alj,Gai_j,N three-layer SLs chosen here. The 
dispersive spectra behave as a series of frequency bands, which are narrower than those in wurtzite 



two-layer GaN-based SLs [20l 134 [35|. The present theoretical scheme and numerical results are 
important and useful for further experimental and theoretical investigations of the polar phonon 
modes behaviors and their effect on optoelectronic properties in wurtzite nitride SLs with complex 
bases. The rest sections of the paper are organized as follows: the phonon states and their dispersive 
equations of a wurtzite GaN-based SL with n-layer complex bases are deduced in section 2 based on 
the transfer-matrix method; the numerical results for the dispersion relation, the electrostatic po- 
tential functions of the phonon modes on a wurtzite GaN/Al^Gai-^N/Al^Gai-yN SL with three- 
layer complex bases are performed and discussed in section 3; Finally, we summarized the main 
results and discussed the significance of the theory described in the current work in the last section. 
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2. Theory 



Let us consider a wurtzite GaN/AL Cl Gai_ Xl N/Al X2 Gai_ a2 N/ .../Al Xll _ 1 Gai_ In l N quasi-2- 
dimensional (Q2D) complex bases SL with n-layer per period (as shown in figure 1). The z axis 
is taken along the [0001] direction (c-axis) of the wurtzite crystals, and the origin of coordinate is 
chosen at the left hand side of some a GaN layer. The widths of GaN layers and the AL^Gai-^N 
layers are di (i = 0, 1, 2, 3, . . . , n — 1), respectively. Thus the period of the complex bases SL is 
L = Yli=o di ■ For convenience, we label the GaN layer with origin of coordinate and its right 
AL-Gai-a-N (i = 1, 2, . . . , n — 1) layers as the unit cell 0. The periodic unit cells of the SL on 
the right hand side of the unit cell are labeled successively as unit cell 1, 2, 3, ... , while those 
on the left hand side of unit cell are labeled as unit cell —1,-2,-3,.... Due to the anisotropy 
of wurtzite crystals, the polar phonon frequencies and the dielectric function becomes direction 
dependent. If one denotes the direction of c-axis as z and its perpendicular direction as t, the 
dielectric tensor can be written as 



et,i(w) 
ei(w) = ( e M (w) 

e,,i(u) 



(1) 



where 



t.i o 

' or 



' J t,T,i 



{i = 0, 1,2, 



(2) 



Here u> z ^, l<j z ,t, <^t,L and ujf.T are the zone center characteristic frequencies of Ai(LO), Ai(TO), 
Ei(LO), and Ei(TO) modes, respectively. The subscripts i = 0, 1, 2, . . . , n— 1 denote the ith-layer 
AL^Gai-a-N materials in any unit cell TV, respectively. 
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Figure 1. Schematic drawing of a GaN-based SL with n-layer complex bases. Here the z-axis is 
taken along the [0001] direction (c-axis). 

In the case of free oscillations (the charge density po(r) — 0) and using the Maxwell equations, 
the electric displacement vector D satisfies the following equation under the cylindrical coordinates, 
i.e. 



V D 



- — ( —) 

P dp {P dp> 



1 d 2 



p 2 d(p 2 



d 2 



0. 



(3) 



It is obvious that equation ^ is a Laplace equation, and its solution $(r) is the electrostatic 
potential of polar optical phonon modes. In fact, via solving the equation ((3]) in the wurtzite SLs, 
one can get the complete polar optical phonon modes. In order to solve equation ([3]) conveniently, 
we define two functions as 



Ci(w) = 
7»M = 

7»M = 



sign [e z< i(u))e t> i(u)] 



\Je t ,i(u)/e z ,i(u), (i = 0, 1, 2,..., n- 1), 

V l e MH/ e »,tHI> CtH ^ o, 
±i v / |^(^y7eI~H|, < o. 



(4) 
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In terms of the theories of two-order differential equations, we know that, when d(uj) < 0, the 
electrostatic potentials of phonon modes in the ?th-layer material correspond to the oscillating 
waves. On the contrary, they correspond to the decaying ones when Q(uj) > [US]. Taking into 
account the oscillating or decaying waves in n-layers of some unit cell in the complex bases SL 
structures, it can be confirmed that 2" types of polar phonon modes are likely to exist in wurtzite 
GaN-based SLs with n-layer complex bases. They are the 10 phonon mode (£i( w ) > 0), the PR 
phonon mode (Ci(w) < 0) [25ll3^. and the (2" — 2) types of the QC phonon modes (at least one 
Ci(oj) > 0, and one Q(lu) < 0, i ^ j) [34, 37, 39| in wurtzite SLs, respectively. These QC modes 
are distinguished and labeled as follows. In some unit cell of the complex bases SL, if the function 
Ci(tjj) takes negative values in j = u, v, w, . . . (j < n) layers, then the QC mode is expressed by 
Q(-m,v,w,... _ p or exam pi ej m a three-layer complex bases SL, the QC modes are denoted as QC 1-11 
if the functions £i( w ) (* = L 2) are negative, and Cs( w ) > 0. 

Next, let us discuss the phonon states and dispersive equations of all the phonon modes in 
wurtzite nitride SLs with n-layer complex bases in a uniform method. Considering the translational 
invariance in p-plane (the perpendicular plane to z axis) of the SL quantum structures, the phonon 
electrostatic potentials in the iV-th unit cell of the wurtzite GaN-based SLs with complex bases 
can be written as follows, 



<$> N {r) = <$> N {p,z 



= e ikt ' p 



(5) 



and 



' <f> N ,o(z), NL<z^NL + d Q , 



4>n(z) = < 



(6) 



[ 4> N ,n-i(z), NL + Y™Jdi<z^(N + l)L 

( A exp[7 fc t (z - NL)] + B cxp[- 7o fc t (z - NL)}, 



e iqNL x < 



Ai exp[-f l k t (z - NL - J2\=a d i)] 
+Bicxp{-jik t (z - NL - YhZo d i)], 

4-i exp[7„_ifc t (z - NL - E"=o 2 d i)] 
k +B„_i exp[-7„_ifc f (z - NL - YJiZa d i)h 



where k f is a wave-vector of phonon modes in the p-plane. And q is a wave vector that will arise 
in the dispersion relation for the collective excitations in SLs. In fact, the phase factor e lqNL is 
introduced to fulfill the Bloch's theorem: 



41|,|42| 



4>(z) = d« z U q {z), 
U q (z + NL) = U q (z), (JV = 0,±1,±2,...). 



(7) 



Via equation (U), it is seen that 7i(w) = ^J\et,i{uj)/ e z ^{ui)\ as Ci( w ) > 0- The potential 4>N,i(z) = 
Ai exp[yiktz] + i3i exp[— 7 ifc t z] in the i-layer of unit cell N behaves as a decaying wave based 
on equation ((5]) in this situation. On the contrary, ji(oj) = ±iy/\et t i(u>) / e z ,i(ui)\ as d(uj) < 0. 
Under this situation, based on equation it is found that the potential function 4>N.i{z) — 
Ai exp[i | 7 i | k t z]+Bi exp[— i | 7 j| k t z] in the i-layer of unit cell N, which just behaves as an oscillating 
wave. This illustrates that, whether Q(u>) is positive or negative, the equation ^ gives the complete 
phonon electrostatic potentials of the phonon modes including the IO modes, the PR modes and 
the QC modes in wurtzite complex bases SLs. 

The phonon potential functions and the normal components of electric displacement should 
be continuous at each heterostructure interface. These are referred to as continuity boundary 
conditions (BCs). Using these continuity BCs, one can obtain a series of BC Eqs. © in GaN-based 
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SLs with complex bases, i.e., 



£z.O 



(f) N ,o(NL) 



4>N-l,n-l(NL) , 

t>N-l,n-l(z) 



NL 



■>N 



dz 



dz 



-NL 



>N 



^ l+1 (NL + J2 d i)- 



')N.i(z) 



dz 



i=0 

z=NL+J2lZ(, d i 
n-2 

i=0 



e z,i+l" 



i=0 

*>N-l,i+l(z) I 

dz \z=NL+J2lZo< 

n-2 



9 AT.- 



£z.n-2" 



>N 



,n-2{z) 



dz 



£-z,n— 1" 



(8) 



Substituting equation ([5]) into equations ((5]), and defining two 2x2 matrixes equations ([5]), (QI 
and a vector (jlip . the equations ([5]) can be expressed in concise forms. The matrixes and the vector 
are defined as follows: 



e Zii exp(7ifc t di) -e z<i exp^ktdi) 



and 



1 1 



(9) 
(10) 

(11) 



Using the above definitions © to (|TTj) . the equations ((HJ can be expressed compactly by the 
following matrix equations (|12p , i.e., 



exp(-igZ)M„_iC n _i = N C , 



M,C 2 = Ni+iCi+i, (t = 0, 1,..., n-2) 



(12) 



(13) 



M n _2C n _2 — N„_iC„_i . 
Via solving inverse matrix on equation (fT2"]l . it is easy to get the relation: 

n-l 

exp(-ig J L)M„_iC„_i = Y[ DjC n _i , 
where 

D, = MTi Nj+1 . (14) 

The condition of the linear and homogeneous equations (|13[) having nonzero solutions for C n -i 
will result in 2 x 2 coefficient determinant equal to 0, i.e., 



Det[exp(-igL)M„_i - JJ Dj] = 0. 



(15) 



This equation (|15[) is just the dispersive equation of phonon modes in wurtzite nitride SLs with 
complex bases. Substituting the dielectric functions and the structural parameters di into the 
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dispersive equation (1151) , the dispersive frequencies of phonon modes in the complex based nitride 
SLs can be worked out as a function of the wave-number k t . It is worth to mention that, as the 
materials become isotropic (i.e. e t = e z ), the dispersive equation (|15[) of anisotropic SL crystals 
will reduce to those of isotropic SLs |43j. In the situation, there are two types of phonon modes, 
i.e., the 10 and the longitudinal-optical- like modes exist in the isotropic SL systems pH I22I |43|. 

3. Numerical results and discussion 

In section 2, we have derived the phonon states and their dispersive equations of all the polar 
optical phonon modes in wurtzite nitride SLs with n-layer complex bases. However, the corre- 
sponding analytical formulas are complicated. In order to get a clear picture for the dispersive 
behaviors of the types of polar optical phonon modes in wurtzite complex bases SL systems, an 
application of the present theories based on a nitride SL with GaN/Alo.isGao.ssN/AIN three-layer 
complex bases is performed, and a numerical calculation of dispersive relationship is carried out. 
The material parameters used in our calculations are listed in table 1 |26l . 1271 . 1.3 ll | . 

Table 1 . Zone-center energies (in meV) of polar optical phonons, dielectric constants of wurtzite 
GaN, Alo.15Gao.85N and A1N materials [2a I27I |33|. 





<*>t,T 










GaN 


69.25 


91.83 


65.91 


90.97 


5.35 


A1N 


83.13 


113.02 


75.72 


110.30 


4.77 


Alo.15Gao.s5N 


71.332 


90.009 


67.382 


93.87 


5.263 



Table 2. Signs of the functions [i = 1 (GaN), 2 (Alo.15Gao.s5N) and 3 (A1N)] and the 

frequency ranges of the phonon modes in wurtzite three-layer GaN/Alo.isGao.ssN/AIN complex 
bases SLs. 
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Modes 


10 


QC 


QC"' 


QC" 


10 


QC" 


10 


QC 


10 


QC" 


10 


QC" 


10 



Before analyzing the dispersive properties of phonon modes, let us first discuss the frequency ran- 
ges of the possible phonon modes in wurtzite SL systems with three-layer GaN/Alo.i5Gao.85N/AlN 
complex bases. In terms of the material parameters given in table 1, the signs of the functions Ci( w ) 
[i = 1 (GaN), 2(Alo.i5Gao.85N), 3(A1N)] and the frequency ranges of the phonon modes in wurtzite 
GaN-based SLs with three-layer complex bases are shown in table 2. From the table, it is observed 
that the characteristic frequencies of GaN, Alo.15Gao.85N and A1N materials divide the frequency 
to axis into thirteen subranges. For convenience, these subranges are labeled as range I, II, . . ., XII, 
XIII from left to right. In each frequency subrange, the signs of the characteristic functions Ci( w ) 
[equation ([2])] are defined. Based on the properties of different types of phonon modes, namely 
behaving as decaying (vibrating) waves in the complex bases material layers of the SLs, the possi- 
ble frequency ranges of these phonon modes in GaN/Alo.isGao.ssN/AIN complex bases SLs can be 
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recognized. It is found that the 10 modes of the system are likely to exist in the subranges I, V, IX, 
XI and XIII due to the positive values of Ci(^) (i = 1, 2, 3) in these five subranges. The PR phonon 
modes do not appear in the GaN/Alo.isGao.ssN/AIN three-layer complex bases SLs because there 
are no frequency subranges satisfying simultaneously the conditions Ci(k>) < (£ = 1, 2, 3). It is 
found that four types of QC phonon modes, namely, the QC 1 , QC 1,U , QC 11 and QC m phonon modes 
are likely to appear in the structure. Here the superscript "i, ii, hi" of QC 1 ' 11 ' 111 modes represent 
GaN, Alo.15Gao.85N and A1N in some unit cell of the nitride SLs, respectively. In fact, the QC 1 
modes appear in ranges II and VIII, the QC 1 ' 11 mode just exists in subrange of III, the QC 11 modes 
also appear in two subranges of IV and X, and the frequency of QC 111 modes falls into the subranges 
VI and XII. Taking into account the condition of solution [44|, the dispersive equation (fT5)) for 
10 phonon modes has no solution in the subranges I, VII and XIII. Thus, the 10 phonon modes 
just appear in the three frequency subranges of V, IX and XI. And the subranges I, VII and XIII 
are termed as the forbidden frequency subranges. Apart from the forbidden frequency subranges 

I, VII and XIII of 10 modes, the dispersive equations of other QC phonon modes have solutions in 
corresponding frequency ranges. Considering the forbidden subrange VII, the frequency subranges 

II, III,. . ., VI are termed as low-frequency range, while the subranges VIII, IX,. . ., XIII are referred 
to as high-frequency range. 



> 




Figure 2. Dispersive spectra of the polar optical phonon modes in low-frequency as functions of 
the free wave-number kt of p-plane in a wurtzite three-layer GaN/Alo.isGao.ssN/AlN complex 
bases SL with d\ — d2 — d-j — 5 nm. The solid lines and dashed lines correspond to the cases of 
qL = and n, respectively. 



Figure 2 depicts the dispersive properties of the polar optical phonon modes in low-frequency 
range as a function of the free wave-number k t of p-plane in a three-layer GaN/Alo.isGao.ssN/AIN 
complex bases SL with d\ = d% = c/3 = 5 nm. The solid lines and dashed lines correspond to the 
cases of phase qL — and ir, respectively. The areas between a solid line and an adjacent dashed 
line are shaded. In fact, these shaded areas are the dispersive frequency bands of corresponding 
phonon modes in SLs. From the figure, it is clearly seen that four characteristic frequencies cj 2i t,2 > 
u>t,T,i ■ ^t,T,2 and w 2j t,3 (dashed and dot lines) divide the low-frequency range into five subranges. 
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In terms of the order of increasing frequency, the QC 1 , the QC'~ n , the QC U , the 10 and the QC 111 
phonon modes appear successively in the five frequency subranges, respectively. These observations 
are completely consistent with the above analysis in table 2. Apart from the 10 phonon modes, 
all the other phonon modes have infinity phonon branches in each frequency subrange. But only 
two branches of 10 phonon modes exist in the low-frequency range, and they are labeled by IOl 
and IOh shown in the figure. For clarity, no more than six branches of QC phonon modes are 
plotted in each frequency subrange of QC modes here. These QC phonon branches are labeled by 
QC 1 / 1 "/"/"^ (i — i ; 2, 3,...), respectively. The subscript i in symbol QC^ 1 denotes the 

order of these QC phonon branches. In figure 1, all the other dispersive curves are monotonous and 
decreasing functions of k t except for the IOl modes and QC" phonon branch. With the increase of 
k t , the frequency bands (shaded areas) of phonon modes become narrower and narrower. And the 
solid lines and corresponding dashed lines converge to the same frequency values as k t d\ approaches 
10. The profound physical reason for this feature lies in the obvious fact that, as the wave-number 
is large enough, the wave-length of phonon is quite short, and the complex bases SLs can be 
looked at as a series of independent three-layer GaNyAlo.15Gao.85N/AlN QWs. Thus, the phonon 
frequency bands of wurtzite SLs reduce to the phonon branches of wurtzite QWs. For example, the 
frequency of IOl modes of the SLs approaches 71.67 meV as k t d\ — > 10, which is just the same as 



the 10 phonon limited frequency in GaN/AIN QWs [25, 27, 29]. Similarly, the QC 1-11 modes in the 



SLs will reduce to the analogous behaviors of PR phonon modes in GaN/Alo.isGao.ssN QWs [26]. 
In a certain frequency subrange of some phonon modes, the higher is the order i of QC^ 1 
(i = 1, 2, 3, . . .), the weaker is the dispersion of the phonon modes. This is a general feature of 
phonon modes in low-dimensional quantum structures [26l 27 . 3~3 |. An interesting feature is that, 



as the frequencies of QCJ~ n phonon frequency bands are lower than the characteristic frequency 
w z,t,2 , these QCJ -11 phonon bands degenerate into the corresponding QC \ phonon frequency bands. 
Also, the QC" 1 degenerating to IOh is clearly observed. But the QC 11 phonon modes are cut off 
at the frequency uj t ^,i , which means that the number of the phonon branches decreases with the 
increase of kt ■ These are the typical features of phonon modes in anisotropic wurtzite confined 



quantum systems |2a427l l34j. The extension of frequency bands (shaded areas) for the QC phonon 
modes is more distinct relative to those for the 10 phonon modes in the low-frequency range. 

The dispersive frequencies uj of the polar optical phonon modes in high-frequency range as 
functions of k t in a three-layer GaN/Aln.isGan.ssN/AIN complex bases SL are plotted in figure 3. 
The meanings of all curves and symbols in the figures are the same as those in figure 2. Same as in 
figure 2, four characteristic frequencies k>t,L,i , w z,l,2 , ^t,\^,i and w Zi l,3 (dashed and dotted lines) 
divide the high-frequency range into five subranges. The QC 1 modes and QC 1 modes cover the 
lowest- and the highest- frequency subranges of the high-frequency range, respectively. The QC 11 
phonon modes occupy the middle subrange of [w z> l,2 , ^,1^2] ■ These are the same as in figure 2. But 
the 10 phonon modes of high-frequency range appear in two subranges, which is quite different from 
those in figure 2. Except for the IOh , the other dispersive curves are monotonous and increased 
functions of k t . As k t d\ — > 10, the 10 phonon modes in the two subranges converge to three 
certain frequency values, which are the same value of 10 phonon modes in GaN/Alo.isGan.ssN/AlN 



QWs [25Ll27Ll29|. The physics of this feature is analyzed above in figure 2. Comparing the frequency 
bands of phonon modes in the five subranges, it is found that the frequency band extensions of the 
10 modes in contrast to those of the QC modes in high-frequency range are quite obvious. This 
is just opposite to the case of the low- frequency range in figure 2, in which the frequency band 
extensions of the QC modes are more distinct than those of 10 phonon modes in low-frequency 
range. In fact, the frequency band extension of phonon modes is due to the periodic crystal structure 
of SLs, which is quite similar to the situation of energy band of periodic crystal structures (45| . 
As the wurtzite GaN/Al x Gai_ x N/AlN SL reduces to the finite-layer GaN/Al^Gai.^N/AlN QW 
structures, the QC (10) phonon bands in wurtzite SLs will naturally reduce to the dispersive curves 



of QC, PR, half-space (HS) and 10 phonon modes in wurtzite QWs [25H2 

As stated in the Introduction, the complex bases SLs possess more adjustable structural pa- 
rameters, and excellent miniband and minigap properties. These characteristics utilized in infrared 
photodetectors, effective-mass filtering, and tuning of the tunneling current 

EMI. 

As the first 
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Figure 3. Dispersion frequencies lj of the polar optical phonon modes in high-frequency range 
as functions of k t in three-layer GaN/Alo.isGao.ssN complex bases SL structures. The meanings 
of all curves and symbols in the figure are the same as those in figure 2. 



step to investigate the polaronic effect on the optic and electronic properties in GaN-based SL 
with complex bases, the phonon modes and their dispersive spectra of the structures are neces- 
sary. As a new and important topic, the geometrical-parameter increases the polaronic effect of 
electronic and optoelectronic characteristics in GaN-based SLs with complex bases and will be 
studied and reported in the future. The 10 and QC modes in binary GaN/AIN SLs were observed 
in experiments [37-39J. The experimental observations on phonon modes in GaN-based SLs with 
complex bases have been rarely reported by now due to the difficulty of crystal growth. With 
the great advance of semiconductor technology, the phonon modes are anticipated to be found in 
the GaN-base SLs with complex bases by the analogous techniques as those used in (37l - l39| . At 
last, we should point out that the present theoretical scheme and numerical results are based on 
non-retardation limit for simplicity. In fact, this treatment is widely adopted to be dealt with the 



crystal dynamics in GaN-based quantum structures [25l - l28l . l33l . |34| . The experimental results of 
angular dispersion of polar phonons and Raman scattering have been proven to be in good agree- 
ment with the calculations based on the non-retarded DCM (38l . l39l |. The recent calculation of 
polaronic binding energy in GaN nanowire also well agrees with relative experimental values [4rj| . 
This shows the reliability and meaningfulness of non-retardation treatment for the phonon modes 
in GaN quantum structures. Of course, a general case of retarded modes can be considered within 
the framework of transverse magnetic (TM) polarization [33] . 



4. Conclusions 

In the present paper, we have analyzed and discussed the completing polar optical phonon 
modes in a wurtzite GaN-based SL with arbitrary-layer complex bases using the transfer-matrix 
method. Via solving the Laplace equations in wurtzite crystal, it is proved that 2" types of polar 
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phonon modes including the PR, 10, QC 1 / 11 / 111 '" phonon modes are likely to exist in wurtzite 
nitride SLs with n-layer complex bases. The analytical phonon states of these phonon modes are 
obtained by means of the DCM and Loudon's uniaxial crystal model. Based on appropriate BCs, 
the concise dispersive equations of these phonon modes in wurtzite complex bases SLs are derived. 
An application of the present theories to a wurtzite SL with three layer GaN/Alo.isGao.ssN/AIN 
complex bases is performed. Numerical calculations on dispersive properties are carried out. Our 
results reveal that there arc five types of phonon modes, i.e., one type of 10 mode and four types of 
QC modes coexisting in the three-layer GaNyAlo.15Gao.85N/AlN complex bases SLs. The dispersive 
spectra of these phonon modes in SLs extend to be a series of frequency bands. It is also pointed 
out that, if the SL structures degenerate into corresponding finite wurtzite QWs, the QC and 10 
phonon bands in wurtzite SLs will naturally reduce to the PR, 10, HS and QC phonon modes 
in wurtzite QWs (2oT - [29l |. The present theoretical and numerical results are important for further 
analysis and discussion of the dispersive spectra of phonon modes and their polaronic effect on the 
optical and electronic properties in wurtzite GaN/Al^Gai-^N SLs. We hope that the present work 
will stimulate further theoretical and experimental investigations of phonon properties, as well as 
their effect on the physical properties in wurtzite GaN-based SL systems. 
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Ct3h 1/1 no/isipHi/ix onTMHHMX cfjOHomB Ta TxHi flucnepciwHi cneKTpi/i 
BK)pi4MT HiTpuflHoV cyneprpaTKM 3i CK/iaflHMMM 6a3i/icaMi/i: MeTOfl 
TpaHCcpep-MaTpi/iLji 

J]. >Ka H r 

ricxniTexHiKa TyaHH>Koy riaHbio, TyaHH>Koy, HapoflHa Pecny6jiiKa Ki/iTaw 

Ha ocHOBi fliejieKTpwHHoT HenepepBHoT MOflejii i MeTOfly Tpa Hccpep-MaTpn L^i , flocnifl>KeHO flonoBHem CTaHi/i 
nojiapHMX onTUHHUx cpoHOHiB y BKDpm/iT GaN-6a30BaHiR cyneprpaTL4i (CT) 3i CK/iaflHWMM 6a3wcaMW. floBe- 
fleHO, mo 2™ Twnw cpohohhwx MOfl icHytoTb y BtopL4HT HiTpuflHii/i CT 3 n-iuapoBWMM CKJiaflHMMW 6a3wca- 
mh. OTpuMaHO a Ha/iiTM h Hi cpoHOHHi CTa h 1/1 L41/IX MOfl Ta TxHi fli/icnepciw Hi piBHflHHfl y CTpyKTypax BK3pi4l/lTy 
GaN/AUGai-^N CT. 3flmcHeHO HucnoBi po3paxyHKM Ha Tpi/iwapoBi/ix GaN/Alo.isGao.ssN/AIN CK/iaflHux 
6a3wcax CT. Pe3y/ibTaTH noKa3yK>Tb, mo e iHTepcpeWcm onTi/iHHi (IO) cpoHOHHi MOflu Ti/ibKW oflHoro Tuny i 
HOTwpw Twnw KBa3io6Me>KeHwx (KO) 4>°h°hhmx MOfl y Tpi/iLuapoBi/ix GaN/Alo.isGao.85 N/AIN CKJiaflHwx 6a- 
3ncax cyneprpaTOK. flwcnepcmm cneKTpi/i 4>ohohhhx MOfl b CK/iaflHwx 6a3wcax cyneprpaTOK po3LuwpioioTbC3 
flo HM3KH HacTOTHMX 30H. CnocTepiraeTbcn , 1140 noBefliHKa KO MOfl peflyKyeTbca flo IO MOfl. U,n TeopeTWHHa 
cxeMa i HwcnoBi pe3y/ibTaTW e flocwTb Kopi/icHHMi/i fliia aHa/ii3y flucnepciwHux cneKTpiB flonoBHeHHX cpohohhwx 
MOfl i TxHix no/iapoHHMX Bnfi \ab\b y BKDpL4HT Ga N-6a30BaHi/ix cyneprpaTKax i3 CK/iaflHUMM 6a3ncaMW. 

K/iiOHOBi c/iosa: omiAHHi cpoHOHHi MOfti/i, AiAcnepcitiHi cneKTpn, cyneprpaTKa Btopuwjy, Merofl 
TpaHcdpep-MaTpwu,! 
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